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Abstract. We introduce a simple and convenient fitting 
formula for the thermal annihilation line from pair plas- 
mas in cosmic sources. The fitting formula is accurate to 
0.04% and is valid at all photon energies and temperatures 
of interest. The commonly used Gaussian line profile is not 
a good approximation for broader lines. 
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1. Introduction 

For two reasons a need has developed to be able to com- 
pute the annihilation line emissivity in a thermal plasma 
both exactly and rapidly. First, relatively broad and re- 
solved annihilation lines have been observed in the galactic 
black hole candidates Nova Muscae (Gilfanov et al. 1993), 
and IE 1740.7-2942 (Bouchet et al. 1991) as well as from 
a source observed with HEAO-1 A-4 (Briggs 1991). With 
some exceptions (e.g. Maciolek-Niedzwiecki & Zdziarski 
1994) a Gaussian line profile is normally used when fit- 
ting the observations. A Gaussian profile does not have the 
correct shape when the temperature becomes sufficiently 
large. Future gamma ray missions such as INTEGRAL will 
be able to determine the annihilation line shape to such ac- 
curacy that the exact annihilation line shape must be used 
when fitting the line. Second, the theoretical modelling of 
radiative transfer in the hot thermal plasmas thought to 
be responsible for the X-ray and 7-ray emission from e.g. 
active galactic nuclei, galactic black hole candidates, and 
gamma-ray bursts requires a fast way to evaluate the ther- 
mal annihilation line emissivity. 

The thermal annihilation line emissivity is exactly de- 
termined by a single, one-parameter integral first derived 
by Svensson (1983, see also Dermer 1984). Svensson (1983) 
also gave some less convenient but very accurate poly- 
nomial fits to that integral. Our aim here is to provide 



a single, simple, but accurate fitting formula valid at all 
temperatures and photon energies of interest. This fitting 
formula can then be used when modelling observations. 
It also allows rapid evaluations of the annihilation line 
emissivity in radiative transfer calculations. 

The structure of the paper is as follows. In § |^, we sum- 
marize previous theoretical results regarding the thermal 
annihilation line emissivity and in § ^ we introduce a very 
accurate fitting formula for the line shape. Finally, in § ^ 
we summarize our results. 

2. Theory of the thermal pair annihilation line 

Consider a plasma consisting of electrons and positrons 
of densities rt_ and n+, respectively. The electrons and 
positrons are assumed to have Maxwell-Boltzmann energy 
distributions at the same temperature, T. Svensson (1983) 
showed using simple detailed balance arguments that the 
emission rate, {dn/dt){x)dx cm~'^ s~^, of annihilation pho- 
tons of dimensionless energy, x = hv/rricC^, in the energy 
interval dx in such a plasma is given by an integral over 
the cross section for the reverse process, photon-photon 
pair production: 

dfi 

— {x,9)dx — n^ri-cdxx 
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Here, K2{l/9) is the modified Bessel function of second 
kind of order 2, the dimensionless temperature 9 is de- 
fined as 9 = kT/nricC^, and a^^[s) is the cross section 



for photon-photon pair production, 
Rohrlich 1976), 
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with Xcra being the photon energy in the 
center-of-momcntum frame. Near pair production thresh- 
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old, s — 1 <IC 1, we have CT^^(s) = <7t(3/8)(s - 1)^2, 
while far above the threshold, s ^ 1, we have crj^{s) = 
c7T(3/8s)(ln4s- 1). 

Dermer (1984) developed a general theory for produc- 
tion spectra from binary collisions in thermal plasmas by 
reducing a six-dimensional integral over the differential 
cross section to a double integral. Dermer applied this re- 
sult to the pair annihilation process, in which case the 
double integral reduces to a single integral over the pair 
annihilation cross section: 



— (x,0)ax = n+ri-cax , ^ 

dt"- ' ^ + e[K2{i/e)f 



rfs2(s-l)aan„(s)e-^/^^ 



(3) 



Here, (Tann(s) is the pair annihilation cross section, and 
s = 7^^, with being the particle Lorentz factors in the 
center-of-momentum frame. We have rewritten Dermer's 
integral in terms of s instead of his chosen integration 
variable, 7r = 2s + 1, i.e. the Lorentz factor of one of the 
particles in the rest frame of the other. The cross sections 
for the annihilation process, (Tann, and its inverse, a^-y are 
related to each other due to microscopic detailed balance: 
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Using this relation we find that the two expressions (1) 
and (3) for the annihilation emissivity, obtained using two 
quite different methods, are identical. 

We now need too evaluate the Bessel function factor 
as a function of and the integral as a function of the 
parameter x0. To an accuracy of 0.06% the Bessel function 
factor can be approximated with (Svensson 1983) 



[1 + 2.00496'-^ -I- 1.47746'-2 + 7r(26»)-^] 



(5) 



which has the correct analytic behavior for ^ 1 and 6 » 
1. Using the two asymptotic forms for the cross section 
at s <C 1 and s ^ 1 allows the integral to be solved 
analytically in the two limits x9 <^ 1 and x6 ^ I giving 



/OO 



aT ^x9 {In AtjxO - 1) ; x9 ^ 1 



(6) 



where rj ~ cxp(— 7e) = 0.56146.., 7e being Eulcr's con- 
stant. Using the asymptotic analytic forms of both the 
Bessel function factor and the integral gives the two fol- 
lowing useful asymptotic expressions for the annihilation 



emissivity 
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At relativistic temperatures, 9^1, the kinetic energy of 

the annihilating pair dominates over the rest mass energy, 
and the annihilation spectrum mimics a Maxwellian with 
a peak at x 9. At nonrelativistic temperatures, ^ ^ 1, 
the rest mass energy dominates and the annihilation line 
peaks at a; ~ 1 + 30/4 with the slight blue shift and broad- 
ening being due to the kinetic energy of the annihilating 
particles (Ramaty & Mcszaros 1981). The nonrelativistic 
expression is valid through out the whole low energy wing 
of the line, and far into the high energy wing up to a pho- 
ton energy of 1/6*. If one only considers the line core near 
x ^ 1, or, more precisely x — 1 <C 1, one gets 



^^^ix,9)dx = n^n.cardx^^, 

a; - 1 «: 1, 



-{x-lf/0. 



«1, 
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which shows that the line core has a Gaussian shape for 
^? ^ 1. It is this nonrelativistic line core approximation 
that is normally used when fitting observed annihilation 
features. Then, however, an exponent of {x — Xcf /9 is 
used with two fitting parameters, the line center energy, 
Xc, and the temperature, 9. 

It is appropriate at this point to briefly mention other 
work on the annihilation line profile. Aharonian, Atoyan, 
& Sunyaev (1983) also derived an single integral expres- 
sion for the thermal annihilation line emissivity, but this 
differs from the self consistent and correct expressions (1) 
and (3) obtained by Svensson (1983) and Dermer (1984). 
While analytical (or numerical) integration of the correct 
expression over photon energy gives the well known an- 
alytical expression for the annihilation rate already ob- 
tained by Weaver (1976), a similar integration over Aha- 
ronian et al.'s expression does not. The reason for the 
discrepancy is the second term I2 in the integrand of their 
equation (A. 4) originating from the term oc cos^ x' their 
equation (33). Using I\ only, gives the correct analytical 
result. 

The multi-dimensional integral for the emissivity that 

Dermer started out with in his derivation has been evalu- 
ated numerically using Monte Carlo methods by Zdziarski 
(1980), Ramaty & Meszaros (1981), and Yahel (1982). The 
results of the first two papers agree with the correct ana- 
lytical results, while the calculations by Yahel give quite a 
different line shape. The reason for the discrepancy is due 
to Yahel confusing the angle between the annihilation pho- 
ton and the particle momentum with the angle between 
the direction of transformation to the lab frame and the 
particle momentum, all in the center-of-momentum frame. 
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The existence of simple analytical expressions for the ther- 
mal annihilation emissivity has now eliminated the need 
for the time-consuming Monte Carlo calculations appear- 
ing in these papers. 

3. An excellent fitting formula 

We now attempt to provide fitting formulas for the inte- 
gral in equation (1) valid for all values of x9. As the anni- 
hilation line is swamped by bremsstrahlung emission for 6 
larger than about 3 (Svensson 1982, Maciolek-Niedzwiecki 
et al. 1995), a useful formula for fitting observable lines 
would need to be valid for, say, x9 < 20. We therefore 
write the integral as 



dssa^.y{s)e 



1/2 



{x6) 



C{x9) 



x6i < 20 , 



(9) 



where C{y) is the correction factor to the asymptotic ex- 
pression (6) for x0 <gC 1. The annihilation emissivity, equa- 
tion (1) then becomes 
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Fig. 1. (a) The correction factor, C{y), as function oi y = x9 
for y < 20. (b) Relative error in per cent using our first Fade 
approximation, equation (11), for C{y). (c) Relative error in 
per cent using our second Fade approximation, equation (12), 
for C{y). 
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any 9, x9 <20. (10) 



The Bessel function factor is given by the approximative 
equation (5) while C{y) was computed to an accuracy of 
10~^ using integration routines from Numerical Recipes 
(Press et al 1986). The resulting C{y) is shown in the top 
panel in Figure |l|. In order to find accurate fitting formu- 
las for C{y), we use two different Pade approximations 
with increasing accuracy. Fitting the numerical results of 
Figure |^ with the Pade approximation 



Ciy) 



l + aiy + a2y 
1 + biy + b2y'^ 



y < 20, 



(11) 



gave the following coefficients: ai — 7.308, 02 — 0.4946, 
bi — 5.002, and 62 = 0.8075. The error as a function of 
y = x9 is shown in the middle panel in Figure 0. The 
maximum error is about 1.4% for x9 < 20 and 0.5% for 
x9 < 10. Similarly, using the expression 



C{y) 



1 + ciy + C2y^ + c^y^ 
l + diy + d2y'^ + dsy^ 



y<20. 



(12) 



we find the coefficients ci = 6.8515487, C2 = 1.4251694, 
C3 = 0.017790149, di = 4.63115589, ^2 = 1.5253007, and 
c?3 = 0.04522338. The error of this fit is shown in the lower 
panel of Figure where the maximum error is seen to be 
0.04 per cent. 

In rare circumstances in numerical simulations one 
may need to know the emissivity for x9 > 20. Applying 



a correction factor, CR (y), to the asymptotic limit of the 
integral, equation (6) for x9 ^ 1, we have 

poo o 

/ dssa^^{s)e-'/'''^ = aT-x9{\n4T]x9 - l)Cu.{x9) ; 
Ji 8 



x9 >20. 

Then the annihilation emissivity becomes 

dTi 3 

— (a:, 9)dx = n^n_C(7Tdx—x{ln4:r]x9 — l)x 
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Fitting the results from numerically evaluating the inte- 
gral gives 



Cniy) 



l-h2.712y-i 
-7800y-* ; 



55.60?;" 



■ + 1039. 8y~^ 
any 9 , y > 20 



(15) 



to an accuracy of 0.1 per cent. 

Our finale recipe for computing the annihilation emis- 
sivity for x9 < 20 is then to use equation (10) together 
with the fit for K2(9), equation (5), and one of the fits for 
C(y), equations ( [Til ) or (|T^. The total error is at most 0.1 
per cent for the latter case. For those rare cases when one 
is interested in the annihilation emissivity for x9 > 20, 
one should use equation (14) together with the fits, equa- 
tions (5) and (15). The total error is at most 0.15 per 
cent. In Figure |2| we show the dimensionless annihilation 
line shapes, dn/dt{x,9)/n^n-caT, using the recipe given 
above for different dimensionless temperatures, 9. The line 
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Fig. 2. Dimensionless thermal pair annihilation line shapes, 
dn/dt(x, 9)/n+n-C<JT, as function of dimensionless photon en- 
ergy, X = hv/meC^, using the analytical fits in the text for 
dimensionless temperatures, 6 = kT/rucf? = 0.01, 0.03, 0.1, 
0.3, 1, 3, and 10 corresponding to temperatures, T, of about 
5 X 10^ K, 1.5 X 10* K, 5 X 10* K, 1.5 x 10^ K, 5 x 10^ K, 
1.5 x 10^° K, and 5 x 10^° K, or to kT of about 5 keV, 15 keV, 
50 keV, 150 keV, 500 keV, 1.5 MeV, and 5 MeV, respectively. 



core is approximately Gaussian for temperatures less than 
about 10^ K, but becomes strongly asymmetric for larger 
temperatures. 



4. Summary and Discussion 

We have derived very convenient formula to use when fit- 
ting observations of thermal annihilation lines or for the 
rapid evaluation of the annihilation line emissivity in ra- 
diative transfer calculations. 

The main result is equation (10) to be used together 
with equations (5), and ( pi] ) or ( p^ ) when computing the 
annihilation emissivity for x6 < 20. Using equation ( p2[ ) 
gives an error of less than 0.1 per cent. The error can be 
reduced to less than 0.04 per cent by evaluating K2{l/9) 
more accurately using, e.g., routines in Numerical Recipes 
(Press et al. 1986). For x6 > 20, one should use equations 
(14) together with equations (5) and (15). The total error 
is at most 0.15 per cent. 

There are a few limitations to the expressions we give. 
They assume a Maxwell-Boltzmann distribution for the 
annihilating particles. The mechanism normally thought 
to be responsible for maintaining such a thermal distri- 
bution is energy exchange through Coulomb scattering. 
This fails for temperatures kT of order nieC^, i.e. 9 > 1 
as discussed by e.g. GhiselHni, Haardt & Fabian (1993). 
Other thermalizing mechanisms may, however, be operat- 
ing. GhiscUini, Guilbert, & Svensson (1988) and Ghisellini 



& Svensson (1990) show that synchrotron self-absorption 
may act as an extremely efficient thermalizing mechanism. 

The expressions given are in the Born-approximation. 
Below temperatures of about 10^ K, Coulomb corrections 
must be included. At such low temperatures, the anni- 
hilation takes place through positronium formation and 
one must also include the positronium continuum due to 
three-photon annihilation (Ore & Powell 1949). 
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